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ABSTRACT 

We report on the optical identification of tfie companion to tlie eclipsing mil- 
lisecond pulsar PSR J1701— 3006B in the globular cluster NGC 6266. A relatively 
bright star with an anomalous red colour and an optical variability (~ 0.2 mag) 
that nicely correlates with the orbital period of the pulsar (~ 0.144 days) has 
been found nearly coincident with the pulsar nominal position. This star is also 
found to lie within the error box position of an X-ray source detected by Chandra 
observations, thus supporting the hypothesis that some interaction is occurring 
between the pulsar wind and the gas streaming off the companion. Although the 
shape of the optical light curve is suggestive of a tidally deformed star which has 
nearly completely filled its Roche lobe, the luminosity (~ 1.9 Lq) and the sur- 
face temperature (~ 6000 K) of the star, deduced from the observed magnitude 
and colours, would imply a stellar radius significantly larger than the Roche lobe 
radius. Possible explanations for this apparent inconsistency are discussed. 

Subject headings: Globular clusters: individual (NGC 6266); stars: evolution - 
pulsars: individual (PSR J1701— 3006B) - binaries: close. 
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1. Introduction 

NGC6266 (M62) is one of tlie most massive and brightest (My = -9.19; Harris 199^) 
Galactic Globular Clusters (GCs), also is characterized by high values of the central density 
(log po ~ 5.47, with po in units of MqPc'^; Beccari et al. 2006). It displays a King-model 
density profile with an extended core (~ 19") and a modest value of the concentration 
parameter (c = 1.5; Beccari et al. 2006). 

Six binary millisecond pulsars (MSPs) have been discovered in M62 (D'Amico et al. 
2001; Jacoby et al. 2002; Possenti et al. 2003, hereafter P03) and it ranks fifth of the 
GCs in wealth of MSPs, after Terzan 5, 47 Tucanae, M15 and M2^. P03 presented phase- 
connected timing solutions, yielding precise celestial coordinates, for three of the MSPs 
in NGC6266. One of them, PSR J1701-3006B (hereafter PSR6266B) displays partial or 
total eclipses of the radio signal at 1.4 GHz near its superior conjunction (in the convection 
adopted throughout this paper, this corresponds to orbital phase = 0.25), clearly due 
to gas streaming off the companion. The pulsar orbit is circular, with a projected semi- 
major axis of only ~ 0.11 i?©. P03 suggested two options for explaining the behaviour of 
the PSR6266B system. The first option is that the pulsar companion is a non-degenerate 
bloated star, whose mass loss is sustained by ablation of its loosely bound surface layers by 
the relativistic wind emitted by the pulsar. In this case, PSR6266B may resemble PSRs 
B1957+20 (Fruchter et al. 1990) and J2051-0827 (Stappers et al. 2001): the optical light 
curve of their companion star presents a maximum when the side of the companion facing the 
pulsar points toward the observer (i.e., at the pulsar inferior conjunction: (f) = 0.75). This is 
a clear signature of the irradiation of the companion surface by the pulsar flux. The second 
option has the pulsar companion as a tidally deformed star overflowing its Roche lobe due to 
the internal nuclear evolution. In this case, the system may be more like PSR J1740— 5340 
(D'Amico et al. 2001; Ferraro et al. 2001), where irradiation effects are negligible (Orosz & 
van Kerkwijk 2003) and the optical light curve of the companion is dominated by ellipsoidal 
variations (Ferraro et al. 2001), i.e. it shows two maxima at quadratures (0 = 0.0 and 0.5). 

We present the optical identification of the companion to PSR6266B, based on high- 
quahty phase-resolved photometry obtained with the Hubble Space Telescope (HST), and 
X-ray emission detected with Chandra . 



^Everywhere in this paper, we refer to the data at jhttp : / / www . physics . mcmaster . ca/ Globular . html 
^See http : //www . naic . edu/ pf reire/GCpsr . html for a catalog of MSPs in GCs with relative references 
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2. Observations and data analysis 

HST OBSERVATIONS. The photometric data presented here consist of a set of high- 
resolution images obtained on August 1 2004 by using the Advanced Camera for Survey 
(ACS), on-board the HST, and retrieved from the ESO/ST-ECF Science Archive. The data 
comprise three i?-band exposures (2 x 340 s and 1 x 30 s exposures), two i?-band exposures 
(of 120 and 340s) and four /7^-band exposures (of 340, 1050, 1125, 1095s). The three 
longest exposures in Ha are the combination of three sub-exposures. An additional set of 
Wide Field Planetary Camera 2 (WFPC2, Prop. 10845, PLFerraro) data has been secured 
through the filter with the specific aim of testing the possible variability of the companion: 

1200 s exposures were taken between May 2 and May 5 2007. In order to best resolve 
stars in the most crowded regions, the Planetary Camera (with a resolution of 0.046"/pixel) 
of the WFPC2 was pointed at the cluster center. The photometric analy sis of the AC S 
dataset has been carried out using the ACS module of the DolphotQ package ( Dolphin 2000 ). 



using the FLT.fits datasets for the photometry, and the drizzled images as references for the 
geometric distortion correction. We set the photometry parameters as recommended in the 
Dolphot manual, obtaining a final catalog of -B, R and Ha magnitudes, calibrated according 
to Sirianni et al. (2005). The final catalog is reported to an absolute astrometric system by 
cross-correlating the stars in common with the data set of Beccari et al. (2006), which has 
been astrometrized by using the standard stars from the new Guide Star Catalog H {GSCII) 
lying in the considered field of view (FoV). At the end of the procedure, the rms residuals 
were of the order of ~ 0.5" both in RA and Dec; we assume this value as representative of 
the astrometric accuracy. 

In order to carefully search for any variability from the MSP companion, we reanalyzed 
the Ha images in a small area around the nominal position of PSR6266B by using RO- 
MAFOT (Buonanno et al. 1983). This package has been specifically developed to perform 
accurate photometry in crowded fields, and it also allows a visual inspection of the quality 
of the PSF-fitting procedure. In particular, we extracted 50 x 50-pixels (~ 2.5" x 2.5") 
sub-images from the ten original Ha ACS frames, each centered on the nominal position of 
PSR6266B. In order to optimize the detection of faint objects, we performed the search pro- 
cedure on the median image. Then, we adopted the resulting mask with the star positions 
to the individual sub-images, and performed the PSF-fitting procedure separately on each 
of them. We adopted the same strategy for the analysis of the nine WFPC2 Ha images. 



■^Only 9 exposures have been used in the analysis since one was made useless by a cosmic ray passage. 
'''The DOLPHOT packages, including the ACS module and all the documentation, can be found at 



http : //pur cell . as . arizona . edu/dolphot/, 
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after the careful alignement with the previously astrometrized ACS frames. The resulting 
instrumental magnitudes were transformed to a common system, and we compiled a final 
catalog with coordinates and magnitudes for all the stars identified in the considered 19 
sub-images. The photometric calibration of the instrumental magnitudes and the absolute 
celestial coordinates were determined by using all the stars in common with the overall ACS 
catalog, calibrated and astrometrized as discussed above. 

Chandra Observations. NGC6266 was observed for 63 ks on May 2002 with the 
Chandra Advanced CCD Imaging Spectrometer (ACIS). The ACIS-S3 chip was pointed at 
the cluster center, its FoV (8' x 8') entirely covering the half-mass radius (1.23'; Harris 1996). 
The data were reduced using the CIAO software (v. 3.3.0), reprocessing the level 1 event files; 
the task acis_run_hotpix was used to generate a new badpixel file. Then we applied the 
corrections for pixels randomization, good time intervals and grades. Our analysis includes 
only photons with ASCA grades of 0, 2, 3, 4 and 6. About 1.5 ks of observations were not 
considered in the analysis because of the high background level. We searched for discrete 
sources in our level 2 event file using Wavdetect (Freeman et al. 2002), performing the source 
detection separetely on the 0.2-1, 1-2, 2-8, and 0.2-8 keV bands. The detection threshold 
was set to 10~^, and the source detection was performed increasing the sequence of wavelets 
from 1 to 16 by factors of \/2. More than 110 sources were detected in the 0.2-8 keV band 
in the entire chip. The astrometry of all the detected sources was corrected by using the 
Aspect Calculator provided by the Chandra X-Ray Center. Then we searched for possible 
optical counterparts both in our HST images and catalog. About 5 bright stars and 2 faint 
and extended objects (likely background galaxies) were found to coincide with the X-ray 
sources well outside the half-mass radius, where the stellar density is relatively low. We 
therefore used these optical counterparts for extending the astrometric solution of HST to 
the Chandra sources, thus obtaining the same accuracy in the absolute astrometry for both 
the optical and the X-ray sources. 



3. Results 

Figure [1] shows the 3" x 3" WFPC2 Ha map centered on the PSR6266B nominal position 
(marked with a cross), as derived from timing (see P03). An accurate photometric analysis 
has been carried out for all the ~ 30 stars found in our catalog within a distance of 1.5" 
from the PSR6266B position. This correspond to 3 times the la uncertainty (0.5") in 
the absolute astrometry of the HST data. In particular, we have extracted 19 individual 
images from our dataset (10 exposures with ACS and 9 with WFPC2) and compared the 
resulting magnitudes for each object in order to search for variability and estimate the typical 
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photometric uncertainties at different magnitude levels. Only one source in the catalog 
(hereafter named COM6266B, marked with a small circle in Fig. [1]) showed Ha variability 
(Ai^Q, ~ 0.2 mag) significantly larger than the typical rms magnitude fluctuations of stars 
of similar luminosity {5Ha ~ 0.02 mag in ACS data, and 5Ha ~ 0.04 mag in WFPC2 data). 
Its celestial coordinates are aj2ooo = 17'^01™12'^.690 and 5j2ooo = — 30°06'48'.'61, whereas the 
mean apparent magnitudes are B = 20.58, V = 19.48, R = 19.02, Ha = 18.65 {V is from 
Beccari et al. 2006). It is located at 0.5" from PSR6266B just at the edge of its positional 
error circle. 

Figure [2] reports the {R, Ha — R) and (R, B — R) Color Magnitude Diagrams (CMDs) 
for all the stars detected in the 20" x 20" box centered on the PSR6266B nominal position. 
COM6266B is marked with a large filled triangle in both the CMDs. As apparent from the 
figure, the star has almost the same luminosity of the cluster Turn Off, but it shows an 
anomalous red color which locates it out of the Main Sequence. The photometric properties 
imply that it is not a degenerate star and that it has a moderate Ha excess, indicating 
the presence of ionized gas surrounding the system. These findings are consistent with the 
irregulaties seen in the radio signal of PSR6266B and supports the scenario that COM6266B 
is a mass losing star. 

In order to firmly assess the physical connection between COM6266B and the binary 
pulsar PSR6266B, we have accurately investigated the time scale of the optical variabilitjo 
(in the 0-1 day range). The Ha time series (19 points) was processed using the GRATIS 
fitting routine, a code developed at the Bologna Astronomical Observatory in order to study 
the periodicity of variable stars (see, e.g., Clementini et al. 2000). Following the procedure 
fully described in Ferraro et al. (2001), the most significant (99% confidence interval)!^ 
periodicity was found at a period of 0.1446 ±0.0015 day, consistent, within the uncertainties, 
with the orbital period of PSR6266B derived from timing (P;, = 0.1445454303(6) day, where 
the figure in parenthesis is the uncertainty, at 99% confidence level, on the last quoted digit; 
P03). Given that, we have folded the 19 magnitude values, by using and the reference 
epoch of the PSR6266B radio ephemeris (Tq = 52047.258199(2); P03). The results are 
shown in Figj3l As apparent, the ACS and the WFPC2 data (obtained about 3 years later) 



^Note that COM6266B has a close companion located at only 4 pixels (^ 0'.'19) West in the ACS map 
(see FigH]). This object turns out to be a normal Sub-Giant Branch star, slightly brighter 0.5 mag in 
B, R and Ha) than COM6266B and not displaying any evidence of variability. Although very close to 
each other, the two stars are clearly separable in the high resolution ACS and WFPC2 images and we have 
verified that the photometric analysis of COM6266B is not affected by the presence of this close star. 

^The significance level has been obtained from the reduced of the Fourier time series fit (at different 
modulation periods) to the Ha data (see for example Figure 2 of Ferraro et al. 2001). 
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agree with each other, drawing a well defined light curve with two distinct maxima at about 
</) = 0.0 and cj) = 0.5. This fact, as well as the phases and the amplitudes of the minima, 
confirms that the optical modulation is associated with the pulsar binary motion and strongly 
suggests that COM6266B is a deformed star overfiowing its Roche lobe (see §1). 

Additional properties of this system can be derived from the analysis of the Chandra 
data. We found an X-ray source located at only ~ 0'.'4 from COM6266B and at 0'.'3 from 
the nominal radio position of PSR6266B. The derived coordinates of the X-ray source are 
Q^j2ooo = 17'^01™12'^.700 and (5j2ooo = — 30°06'49'.'08, and the la uncertainty area on its posi- 
tion is encircled in Fig{T] (dashed line). Given the number of sources (50) detected within 
the half-mass radius of NGC6266, the probability of a chance superposition of a X-ray 
source within 0'.'3 from the radio position of PSR6266B is ~ 7 x 10~^, strongly support- 
ing the association of the detected source with the PSR6266B system. The (background 
substracted) photons counts are IS.StH (soft: 0.2-1 keV), 23.7+4:9 (medium: 1-2 keV) 
and 9.8^3 1 (hard: 2-8 keV), implying HRl^ 1.3 and HR2~ 0.4 (where HRl=medium/soft 
counts and HR2=hard/medium counts). Hence the source is harder than the typical MSP 
population observed in GCs (e.g. Grindlay et al. 2002) and the HRl, HR2 values re- 
semble those of a source in which the pulsar wind shocks the material released from the 
companion (see, e.g., the cases of PSR J1740-5340 in NGC6397, and PSR J0024-7204W 
in 47Tuc; Grindlay et al. 2002, and Bogdanov et al. 2005, respectively). Even if the small 
number of photons prevents a detailed spectral analysis, for an absorption column density 
Nh = 2.2 X 10^^ cm~2 (Pooley et al. 2003), the counts in each band and the HRl and HR2 
values are consistent^ with a power law model having a photon index ~ 2.5, and a total 
unabsorbed fiux Fx ~ 10"^'' erg cm~^s~^ in the 0.2-8 keV band, translating (for a distance 
of 6.6 kpc; Beccari et al. 2006) into a X-ray luminosity of Lx ~ 6 x 10^^ erg s~^. 



4. Discussion 

Many observed features of COM 6266B (namely, its anomalous red colour, the excess, 
the shape of the light curve, the X-ray emission from the binary) are suggestive of a tidally 
deformed star which is experiencing heavy mass loss, similar to the system found in NGC 6397 
by Ferraro et al. (2001). On the other hand, other photometric properties of the source do 
not fit this picture, as described below. 

By assuming the reddening E{B — V) ~ 0.47 and the distance quoted by Beccari et al. 
(2006), and adopting [Fe/H] = —1.1 and an age of 12 Gyr for NGC 6266, we have determined 



^The PIMMS tool has been used for these estimates: 



http : / /heasarc . gsf c . nasa . gov/Tools/w3pimms . html 
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the physical parameters of COM6266B from the comparison of its position on the CMDs 
(Fig. [2]) with the isochrones of Pietrinferni et al. (2004). The resuhing effective temperature, 
bolometric luminosity and radius of this star are Teg ~ 6000 it 500 K, Lboi ~ 1-9 ± 0.2 Lq, 
and i? ~ 1.2 ± 0.2 i?Q, respectively, where the quoted uncertainties are conservative esti- 
mates, essentially due to the uncertanties in the distance modulus (±0.15 mag) and in the 
reddening (±0.05 mag). If the luminosity variations shown in Fig. [3] are mainly due to tidal 
deformations of the companion star, the stellar radius R is expected to be of the order of 
the Roche lobe radius Rl. The latter can be estimated from the mass function and the 
projected semi- major axis of PSR6266B (derived from pulsar timing; P03), provided that 
a pulsar mass Mpsr and an orbital inclination i are assumed. As already discussed in P03, 
the presence of radio eclipses indicates that i is not small. Conservatively taking i ^ 20° 
and assuming M^st = 1.4 Mq, the companion star is expected to span a mass interval rang- 
ing from 0.125 {i = 90°) to 0.41 M© {i = 20°)@, corresponding to Rl ~ 0.26 i?© and 
0.40 i?0, respectively. The value of Rl is only slightly affected by the choice of Mpsr, be- 
ing, for instance, Rl ^ 0.45 Rq for 0.7 M0 < Mpsr < 3.0 M© (a safely large neutron star 
mass range; see, e.g., Freire et al. 2008). Note that in order to match the estimated Rl 
with the observed R it should be i ~ 3° (for Mpsr = 1-4 Mq), thus implying an unreliable 
companion mass of ~ 6M0. In summary, for all plausible binary parameters, it turns out 
that the observed stellar radius of COM6266B is significantly larger [R ^ 3 Rl) than the 
expected Roche lobe radius of the companion to PSR6266B, i.e. COM6266B appears ^ 10 
times brighter than the expected luminosity of any pulsar companion having an effective 
temperature Tg = 6000 K. 

Under the hypothesis that COM6266B is the companion to PSR6266B, what is the 
origin of this remarkable discrepancy? A number of possibilities are examined below. 

(i) The binary system could not belong to the cluster. In this case the adopted distance 
would be inappropriate and the value of the radius deduced from the luminosity would be 
largely overestimated. On the other hand, the celestial position of PSR6266B is at less 
than 2" from the cluster center, and the pulsar shows a negative value of the spin period 
derivative, as expected if PSR6266B lies in the cluster potential well. 

(ii) The optical luminosity of COM6266B could be dominated by strong non-thermal 
processes, possibly triggered by the pulsar spin down power Lsd ~ 10 L0 (see P03). Also 
in this case the value of the stellar radius deduced from the observed luminosity could 
be significantly overestimated. However, At the orbital separation of 1.4i?0, the energy 
captured by a Roche lobe filling companion would be of the order of L ~ 0.1 — O.2L0 which 



Note that in both cases the orbital separation of the system turns out to be a '--^ 1.4 i?Q. 
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is too low to explain the observed value. Moreover, it is hard to explain how a large non- 
thermal component may generate a light curve with maxima that fall at orbital phases nicely 
consistent with those expected for light curves modulated by ellipsoidal variations. 

(in) COM6266B could be the blend of two low- luminosity stars: one star would be 
the optically variable companion to PSR 6266B and the other one a non- variable star in the 
foreground (background). While stellar blends are not uncommon in the highly crowded 
GC cores, this possibility only partially alleviates the problem. In fact, by assuming that 
COM6266B is the blend of two stars of equal mean Ha magnitude (i.e., = 19.35), 
the observed ellipsoidal modulation could only be produced if the pulsar companion varies 
by AHa ~ 0.35 mag, which is a rather extreme value for variabilities induced by tidal 
distortions. Moreover, even in this case, the star would have a radius R ^ 1.7 R^. Matching 
R with Rl would require that COM6266B is e.g. the blend of a non- variable star with 
Ha — 18.7 and a star with mean magnitude H^ ~ 21 varying by AHa > 3 mag. We note 
that similar considerations would apply also if PSR 6266B were a triple system. 

In conclusion, while option (Hi) seems to be a quite unlike ad hoc scenario, options (i) 
and (a) will be finally addressed as soon as spectroscopic data are available for the system. 
This will allow us to investigate the nature (thermal or not thermal) of the optical emission 
and the association of the source to the cluster (via the determination of the radial velocity), 
hence shedding some light on the nature of this puzzling object. 

This research was supported by the Agenzia Spaziale Italiana (ASI 1/088/06/0), PRIN- 
INAF2006 and by the Ministero dell'Istruzione, dcll'Universita e della Ricerca, and it is part 
of the Progetti Strategici di Ateneo 2006 granted by University of Bologna. AP and NDA 
also acknowledge financial support from PRIN-MIUR2005. We thank the Referee for the 
careful reading of the manuscript. 
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Fig. 1.— Finding chart for COM J1701-3006B showing the median-combined WFPC2 
image. The region covers 3" x 3"; the cross marks the nominal position of PSR6266B, with 
the la error circle (0.5", solid line). The smaller circle shows the optical counterpart, whereas 
the dashed line represents the la circle (0.7") in the Chandra position. North is up and East 
is to the left. 
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Fig. 2.— {R, - R) and (it!, B-R) CMDs for stars identified in a region of 20" x 
20" centered at the nominal position of PSR6266B. The optical counterpart to the pulsar 
companion COM 6266B is marked with a large filled triangle. 
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Fig. 3. — Ha light curve of COM 6266B obtained by using the period and the reference 
epoch of the radio ephemeris of PSR6266B. Asterisks represent ACS archive observations 
performed in August 2004; large empty triangles are the WFPC2 data collected in May 2007. 
The phases of quadrature and conjunction of COM 6266B are reported for clarity. The solid 
line represents the Fourier time series (trucated to the second harmonic) that best fits the 
data. 



